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1 Introduction 



Calculation of the light-shifts of both the ground and excited states of the 
atoms are mostly easily performed in the formalism of the so-called "dressed" 
atom (see e.g. [|l|). Within the framework of this approach the shift, A£„, of 
the level n occurring due to the incident photon beam of the energy, u;, and 
polarization, e, is defined as follows: 

P (n\D • e*\r)(r\D • e\n) (n\D • e\r){r\D • e*\n)\ 

/\en\Oj) = 5 m. y < 1 > 

2eon Sc ^ [ En + u; - Er En - OJ - Er J 

(1) 

Here D = er is the dipole momentum operator of the electron; P = ^eqScE^ 
is the power of the photon beam with the cross section S and the electric 
strength E; 3? designates the real part. According to (P, the calculation of 
the shift, AEn, reduces to that of the tensor of dynamical polarizability (DP), 
al^{u), of the state n. It is defined as (see e.g. 0) 

^i7(^,) _ ^j^i HDi\i"){r\Dj\n) ^ {n\Dj\r){r\Di\n) \ 

Jy. I En ~\~ — En — ^ — J 

Here Di denotes the iih component of the vector of the dipole moment. An ah 
initio exact analytical calculation of DP is evidently possible for hydrogenic 
systems only and the corresponding results are known (in principle) for all 
states already for long time (see |^ and references herein). They are usually 
performed by means of the exact explicit expression for the non-relativistic 
Green's function either in coordinate or momentum spaces. As a result, DP 
are expressed in terms of the special (Appel) functions whose complicated 
mathematical structure makes accurate numerical (and analytical) analysis 
of a^^ci;) hard to carry out. Especially this argument refers to the case when 
the photon energy, uo, lies in the vicinity of the threshold: uj ^ In-, In 
being the ionization potential of the state n. It is not surprising, therefore, 
that calculations of this type seem to be available for the hydrogenic Is- 
state only (see |Q). In fact, provided that the principal quantum number, 
77,, of the level is fixed, a^(u;) has singularities when the photon energy, 
uj, is in the resonance with the higher/lower discrete levels of the atom Q: 

^This singularities can be avoided if the finite widths of aU atomic levels are taken into 
account. 
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uo = \ei — en\-, i > n,i < n. Beyond that, DP acquires also a non-zero 
imaginary part if u > In- For such u the real part of DP, ^al^^u), describes 
(as for u < /„) the shift of the level, whereas imaginary part, Q'a^ (to'), allows 
for a decay probability (photoionization) of the atom under the action of 
the photon field. Namely, according to the optical theorem [0: a^'^^uj) = 
47rau;S5Q;^(a;), where ctI^^lu) is the total photoionization cross section of the 
state with the principal quantum number, n, and a = e^/ific) = 1/137 is 
the fine structure constant. In mathematical terms, has unremovable 

singularity at the point u = In, so that a great care must be taken to make 
numerical calculation within this region of u stable and highly accurate. 
The most efficient way of achieving that consists in combining numerical 
methods together with analytical ones. It is relevant to point out that the 
detailed description of the DP's behavior of the mentioned type proves to be 
of particular importance, e.g. for the problem of the light-shifts' calculation 
in muonium atom, denoted (/i"*" — e~)°. This is due to the fact that in the 
highly accurate experimental measurements of the Is — 2s energy splitting 
in this exotic system, which are in progress now, the energies, uJi and 002, 
of two incident photon beams are supposed to be in the resonance with the 
following transitions: Is + hui + hui 2s and 2s + huj2 —* ep (see for 
more details). Therefore precise calculation of the corresponding light-shifts 
of Is- and 2s-levels would be rather desirable. 

In the current letter we present the results of analytical and numeri- 
cal calculation of <y]^{uj) together with the corresponding light-shifts in the 
ns-states, n = 1,2, of the muonium atom. It should be pointed out that 
the similar calculation by Beausoleil employing a pure numerical scheme 
proves to be incomplete. Besides, in the contrast to the usual technique (i.e. 
by means of the Green's function) the current calculation is carried out in 
the fashion of Sternheimer ^ where (exact and analytical) summation over 
the intermediate states r in (|^) is reduced to solution of a certain differential 
equation. Such an approach, which is applied to the problem under consid- 
eration for the first time, to our knowledge, seems to be rather instructive. 
Apart from its self-contained academic interest, it may also give certain ad- 
vantages in treating the higher ns-states (n > 4) of the hydrogenic systems, 
as well as for exact calculation of the various tu-dependent sums of the form: 

{n\\r\\s){s\\r^'\\n) ^ (n| |r| |s) (s| |r^| |n) 1 
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Here /i is an arbitrary number, being not necessarily positive and integer; 
(s||r^||n) denotes reduced matrix element. Such type of expressions make 
their appearance in numerous problems of atomic physics. 



2 Calculation of the light-shifts 

2.1 General consideration 

For the ns-states under consideration Eq. (|I|) can be reduced to the following 
angular- and spin-independent form |^ 

P P 
^ens{uj) = - 2g • e*)5Raf,(^) = - ^ ^aU^). (4) 
2eoh be zEoh be 

Here Qnsi^) denotes the so-called scalar DP (henceforth the atomic units, 
e"^ = h = nie = I are used): 

Sr \ I ^ j {ns\\r\\kp) {kp\\r\\ns) {ns\\r\\kp){kp\\r\\ns)\ , , 

«ns(^) = ~oXl \ T^Ti ^ f' 

^ I -T UJ — Skp -r ^U £ns — UJ — Ekp ) 

which involves the radial integrals only. Summation is performed here over 
complete set of discrete and continuum p-states of the Coulomb field with 
the charge Z. The infinitesimal positive imaginary constant added in the 
denominator of the first term in the sum defines the sign of 9Q;f^(co') occurring 
if u > \ens\, \£ns\ = Z'^/{2n'^) being the ionization potential of the ns-state. 
Calculation of a^^^u), Eq.(^, is actually the final aim of our consideration. 
Let us introduce auxiliary function ^„(r; E) by the equation: 

i.„(r;E)^tMMr}M^ (6) 

E = Ens ± UJ + iO being a parameter. In terms of '?/'„(r; E) DP, a^si^j), is 
expressed as: 

"fs(^) = [{ns\\r\\4jn{£ns + + iO)) + {ns\\r\\iljn{£ns - uj))] . (7) 
By acting on ^ipnir) of Eq.(|^) with the operator, 

2r^ dr i dr ) r ' 
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and by virtue of completeness of the set \kp) one immediately obtains the 
following inhomogeneous differential equation obeyed by the function ipnif)'- 



-2Tr { ^'^1 + 2 - 1 + ^) = 2ri?..(r). (8) 
r'^ or \ or J \ r / 

Here -R„s(r) denotes the radial non-relativistic coulomb s-function. In terms 
of the new parameters, 

^ 2Z 

P = —r, (9) 



Eq.(H) takes the form: 

?| [p'^) + H - ;^ + ^) = i (i) • 

Finally, on introducing the new auxiliary function ({p) as 

Mr) = pe-'/'ap), (10) 
we result in the equation of the hypergeometric type, 

pC'ip) + (4 - P)C'(P) + - 2)C(P) = \ [^)\e^"Rns (^) . (11) 

Its solution, C(p); is supposed to be subject for 3?z/ > —n {n is fixed) to 
conditions: 



C(p)=0(l), asp^O, C(p) = o 



n + V 
exp ( ^p^ 



as p — s> oo. (12) 



These follow directly from the definitions, Eqs.(p|),([T(]|), since (see [0): Rnsij") " 
^-Zr/n^ as r — >• oo, and Rnsij) = 0(1), Rkpij) = 0(r), as r — > 0. It should be 
emphasized that Eqs.(|T^ are consistent with the "orthogonality-condition", 

{2p\ipi} = — , if n = 1, or {np\^pr^} = — , if n > 2, (13) 

rj — -C/ — Enp 

which follows from the Eq.(|^) and by virtue of the orthogonality of the kp- 
functions: 

{kp\\r\\k'p) = Sk,k'- 
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Relations (|T3|) may be used as an additional check of correctness of the func- 
tion ipn{r). 

The general solution of ( [n| ) has the form: 



C(p) = D^^2 - V, 4; p) + ^2^(2 - v, 4; p) + Co(p). 



(14) 



Here $(2 — i/, 4;p), \E'(2 — z^, 4;p) are the regular and irregular solutions 
of homogeneous hypergeometric equation []I0|; Di^D^ are some arbitrary 
constants which will finally be chosen to comply with (0); Co(p) is some 
particular solution of Eg. (pA]) . For the distinguished case: v = 2, explicit 
general solution of Eq. (|TTD reads: 



C.=2(P) 



6p 



{2eP + pe" + p'eP + p'Ei(-p)) + + Co(p). 



Here Ei(— p) stands for the integral exponential function. 
Let us seek Co(p) in the form: 



Co(p) = 7^ / e^'mdt. 



(15) 



Here the integral is taken along some contour 7 in the complex plane of t. 
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Figure 1: The contours of integration in Eqs.(20),(21) and (24). 

It has to be chosen to comply both with (^2|) and the form of the contour 
7i (see below) in the integral representation of the right-hand side of ([TT|) . 
Besides, after passing along 7 the integrand of ( [T3| ) should return back to its 
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initial value. To find ^{t) and establish 7 we can write first (see e.g. 



4 \Z 



Vf>_ 

2Z 



1 
2 



n — v . 

I pexp ( ^^p)<l)(-r2 + 1,2; z/p/n) 

'nZ 



2n 



Then, we use the well-known integral representation [IC] of the confluent 
hypergeometric function, $(— n + 1, 2; vpjn): 



$(-72 + 1, 2; vpin) 



— I exp (yptlnM-ty^il - tYdt, (16) 

mn J-yr, 



2mn J 70 

where the contour 70, shown in the Figure 0(a), is passed in the counterclock- 
wise sense along the circle of an arbitrary radius, e. By means of Eq.(0) and 
after two variable changes the right-hand side of ([TTD can be finally express 
as 



4 \Z) ^ \2Zj 4m 



V 



'nZ 



71 



n — V 



2n 



-n-l 



t - 



n + u 

2n 
(17) 

It has thereby the form similar to Eq.(|15[). Here the contour 71, drawn in 
the Figure |l|(b), is passed in the counterclockwise sense around the point 
{n — u)/ {2n) along a circle of an arbitrary radius. A substitution of ([T5|) and 
(|ll in Eq.(0) yields: 

3 



n-l 



dt. 



t(i-t)r(t)+(2t+z.-i)e(t) = M-^ 

^ V VnZ 



n — V 
2n 



-n-l / ^ _^ z/\"-i 
2rr 



The general solution of this equation can be written as 



'nZ , 



l-t) 



00 

l+v I 



(t- ^] 

-u-2 V 2n J 



(18) 

Here Cq is the arbitrary constant which can be set to for convenience. The 
integral is assumed to be taken along any path in the complex t-plane which 
does not pass through the points: {0, 1, (n — z/)/(2n)}. 

If combined with (|T3p, Eq. ([TB|) defines desired particular solution, Co(p)) 
of Eq.(0): 



dt. 



Co(p) 



1 1 



2m 2 



'nZ , 



l+u 



00 1 1 



\l-t)- 



2n 



n-l 



n—v 
2n 



n+1 

(19) 



dt. 
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In view of the given above argument, the contour 7 in (]T9|), is still a free 
"parameter" , provided that it is deformable into 71 . Let us show how we can 
fix it. In fact, the integrand in Eq.(^^ is analytical function in the t-plane 
cut along any path with the ends at t = 00, = {n — z/)/(2n), the latter 
being its logarithmic branching point and the pole of the nth order at the 
same time. Owing to that, we can choose 7 to be a curve which starts at —00 
at the lower edge of the cut, runs along the real axis, encircles the point to in 
the counterclockwise sense and runs back to —00 along the upper edge of the 
cut (see Figure |](c)). Such a contour is topologically equivalent to 71 in the 
Figure |l](b). Moreover, after passing along 7 the integrand returns back to 
its initial value, since it decreases exponentially as t ^ —00. Here we have 
temporary assumed that 3?(i/) > n, '^u = 0. This restriction will however be 
released later by means of the analytical continuation in u. For the contour 
7 under consideration the integral ( [T^ ) can be split into two (independent) 
parts: («) along two edges of the cut and (m) along the circle centered at to- 
The part (i) is reduced in its turn to the integral of the jump at the cut of 
the integrand of Eq.(|T9[). It has the form: 



...dt 



cut 




e-^^(-x)^-^(l + x)^+^x 



X < 



^ n+u 

u-2 \ 2n 



r-2(l-t) 



n-1 



n+1 



dt 



-(...) 




f-'il-t) 



n + u 
2n 

roo 



t=—x — iO 
n-1 



dx 



t=-x+iO, 



t=to 



roo 

X / e-''"(-a;)^-"(l + xf^" dx.{2{)) 

In 

Conversely, the part (m) is expressed as a residue of the integrand at t = to: 




, 00 

1 + !/ / ^v-2i 



ePH^-\\-t) 



Combining Eqs.(|2D|) and ( pij ) together we get: 

3 



0-^ 



n-1 



n+1 



dt 



.(21) 



Co(p) 



1 / V 



1 d'' 



2 V \raZ / n\ df 



X 



i=io 
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X / e 



+ res 

t = tQ 



oo 

l+v I 



-v-2 V 2n 



J- n— 



n+1 



^ .f22l 



By virtue of the integral representation [|T^ of the function \I/(q;, 7; z), enter- 
ing Eq.(|4D, 



1 

^(a, 7; z) = —— / e-"*r-Vl + tV^'^-^dt, > 0, (23) 
r(a) JO 

we can finally write: 

Mr) = ^^^^(2 - z/, 4; p) + D2^{2 - u, 4; p) + Co(p)) • (24) 

Here we have redefined the constant D2 without changing its notation; Co(p) 



is obtained from Co(p) of Eq. (22) by means of the substitution: 

ly — n 

. . .dt ^ — / ^" . . .dt. 



2n 



(25) 



It it clear now that for u : 2 — u ^ 0, —1, —2, . . . the function tpn{r) of 
Eqs.(g2l),(g5|) will satisfy conditions (|g) if we set: Di = D2 = 0, since [0 



$(2-z/,4;p) 
^(2-z/,4;p) 



r(2-z.) 



as p ^ 00 
as p — s> 



P' 



u-2 



p 



as p — s> 00 
^■^5 as p ^ 



We can adopt the same choice of Di also for u = 2,3,4,..., being of no 
physical interest. Whence, one finally obtains: 





V 




nZ 


u — rt 




2n 

e 






res 




t=to 



pe 



[ n! dt" 

^-'{l + xY+'^dx - 



-u-2 



n + u 
2n 



n-l 



t=to 



ft- ^ 



n-l 



n+1 



(it 



.(26) 
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In the given derivation we assumed that parameter v is subject to condition: 
n < 3?(z/) < 2, = 0. Analytical continuation of (^) on all u is achieved 
by means of the substitution: 



i^ — n f — n 



I e~P''{-xY-''{l+xY^''dx fl ' / e-P'^i-xf-^il+xf'^^dx. 

Jo ^ ^ Jo+ 

(27) 

Depending on the cases: 3?z/ > n or 3?z/ < n, the integral in the right-hand 
side here is taken along the paths shown in the Figure 0(a),(b). Each of 
them starts at the point {—to) = {u — n)/ {2n) lying on the lower edge of the 
corresponding cut, encircles the origin in the clockwise (counterclockwise) 
sense, and ends up at (—to) lying on the upper edge of the same cut. In the 
following we shall use for simplicity a sign of the ordinary integral but imply, 
whenever necessary, the substitution (^). Alternatively, the same analytical 
continuation can be achieved by means of the identity: 



e-^"(-a;)^-'^(l + a;)i+"da; 







, . , $1 2-z/,-l-z/,3-z/,^— p , 3?z/<2.(28 
1 — V \ 2n J \ 2n 2n / 

Here <I>i(. ..) denotes degenerate hypergeometric function of two variables 
defined by the following series |jTl|: 



$i(a,/3,7,x,y)= f f N<1. (29) 

^0^0 {l)m+nm\n\ 

In particular cases, n = 1, 2, 3, explicit expressions for the functions ijjnir) 
are combined below. 

n = 1 

^,(r) = -—-l^;^ f^Vpe-^/^ /'^e-^V-^(l +x)i+'^rfx + 



Z3/2(zy2 _ 1)3 ^ zy + 1 



n = 2 
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512^2//^ /Z/-2 



-p/2 



v-2 

4 



2^3/2(1.2 _ 4): 



-i/p/4 



z/(z/2 - 4)p + 4(z/2 + 4) (31) 



V^3(r) 



n = 3 
, + 3 / JO 



n/3i.3 



27Z3/2(z/2 -9)2^'' 



-;/p/6 



z/2(z/^ - 9)p2 - 6u{iy' - 27)p - (306//^ + 486) (32) 



It can be proved (see Appendix A) that Eqs. (|30D -(|32D satisfy orthogonality 
condition of Eq. (|T^ . 

Let us also give for reference explicit expression for the leading term in 
the expansion of ipni'f') in the small parameter u/n ^ 1: 



Mr) X -2 

n>i/ \ \/nZ , 



1/2 _^,fl 

e 

V2 



16z/ / e-P' \--t 



1 

- + 1 ) dt + l 



< 2.(33) 



This identity follows directly from Eq.(P6| 



2.2 Calculation of the light-shifts in particular cases 

Now we are in a position of calculating the matrix elements entering Eq.(0) 
for n = 1,2, being of particular interest. Integration, using Eqs.(P(]|-PT|), 
yields the following equivalent forms of (ls| |r| l^/^i), (2s| |r| |?/'2) (see Appendix 
B): 



;(ls||r||?/'i) 



512Z/9 



Z4(z/2 



,1-u 



+ (-(Si/') 



2z/2 (2z/2 



(z/2 - IJ 



512z/^ 



2F1 4, 2 - z/, 3 - u; 



+ 



V - 1 

Z/+1 
2z/2 ^- -2 



+ 



Z^ iv^ - 1 
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1 2^^u^ /"i t^^" 



1/2 (z/4 - 64z/2 + 112 

220^9 ^1 ^2-. 



^ / /'^-2\2\ 
■ 2^1 4, 2 - z/, 3 - z/; + 



z/2(z/4 - 64z/2 + 112 
^ ' 2i^i 4,3-z/,4-z/;^ + 



Z4(z/ + 2)i2(i/-2)(z/-3) ' ' V ' ' 'Vz/ + 2 

i/2(z/3 + 38z/2 + 84i/ + 
+ 16- 



Z\v - 2){v + 2) 

Eqs.(p4D,(p6D prove to be convenient for an analysis of the DP's behavior 
when uj ^ Z'^j(2r?^ = In and a; ~ 0, as well as for calculation of the DP's 
imaginary part; Eq.(^) is suitable for numerical calculation of the light-shift 
of the 2s- level; Eqs. (|35|) , (|38|) , and (|39D define the matrix elements in terms 
of the hypergeometric function ||T0[, 2-^1 (• • •)• ^^^^ form they easily admit 
analytical continuation in u. 

By means of Eqs.(0),(|) we can explicitly express now DP, afg(co'), Q;f^(c<j), 

as: 

«fsH = -^((i^lkll^i)L=.n+ (i^ll^ll^i)l.=J' (40) 

«fsM = -^((25|kll^2)L=,,,+ (2s||r||V^2)L=J. (41) 

Here we have introduced the following notations: 

Z Z 

T^ll = I 1^12 = I = (42) 

I -2{eis + oj + iQ) J-2{eu-oj) 
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1^21 = ; 1^22 = , =, (43) 

-2{e2s + + iO) J-2{e2s-uj) 



so that the following identities hold true: 



/I 4 r/l 4 ' ' 



Here eis = —Z'^/2, = —Z'^/S are the energies of the Is- and 2s-levels. 
We assume, according to the standard rule of analytical continuation of the 
square root (see 0), that 

' > if E < 



-2(E + zO) \ ^ ifE>0 

Hence, by choosing the matrix elements |r| I"?/;!), (2s| |r| |'?/'2) in the forms 
of Eqs.(BH),(|5S|), DP can be expressed as: 



1 A 512Z/L 



E 



X2Fi(4,2-z/i™,3-i/i^;(^i^^) ),(45) 



7^i2Fif4,2-^2™,3-z.2^;(^^^) ) .(46) 



Eq. (|i5| ) is in agreement with the well known result of Gavrila P , as it should. 
Eqs.(^5])-(^) are rather inconvenient, however, for numerical calculation of 
DP for energies lying above the threshold of the levels, i.e. when param- 
eters i^ii, z/21 become purely imaginary (cf Eq. (^4]) . Besides, they are also 
unsuitable for obtaining various asymptotics of these quantities. As was al- 
ready mentioned above, for these purposes the integral forms of the matrix 
elements, Eqs.(^4|),(^6|),(|37D, prove to be more convenient. Below we com- 
bine various most important results of such calculations (see Appendix for 
details). 



1. The case: cj/Z^ < 1. 

9 1 , 319 ^ a; 

2^4 



^Uu) X + +••■ (47) 
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2. The case: a;/^^ > 1. 

^ , ^ 1 4v/2, , AttZ^ 

w/z2>i cu^ 6u!^ 3 o;"/^ 3 lu^ 

- ^ (-336i - 35i7r^ + 336 + 327r^ - 3(1 + i)7rln(8ZVa;)) 



1 1 V2,_ Z^ .71 Z 
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2 — - 



a)/Z2>l g ^4 g ^9/2 g ^5 

1 

- 2304^ ~ ^1^^' + + + + i)7rln(2ZVu;)) + 

3. Calculation of (a;), a; > /„ = Zy{2n'^). 

647r g~4»;arctan(l/»y) 

W-^) = ^ (1 + ^2)3 .^-I^nl. (51) 

Accordingly, tlie photoionization cross sections, a^J{uj) — 4:Trauj'^a^{(jj) , 
are defined as: 



a = pi r rj = J > 0, (53) 



4. The case: a;/4 ^ 1 + 0. 

X ^-|^e-^Ei(4)(l-i)- 
3 a)//is^i+o 6Z^ 6Z^ 



I/O M M , ^ 18896 131072 o^.,,,,, 

-2srK/.2 ^e-«Ei8 1-^ 

3 3Z* 3Z* 



^ (1016 - 2816e-^Ei(4)(l - t)) \ + . . 



f 709952 - 4849664e-^Ei(8)(l - i)) \ + . 

9Z** ^ '' T]^ 
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Here parameters i] are defined by Eqs.(p3D,(|5^; e = 2.71828 . . .; Ei(. . .) 
stands for tlie integral exponential function; condition 77 1 is assumed 
in either case. 

5. The case: Z^/(2ujn^) < 1; Z,uj are fixed. 

By means of Eq. (|33D and in view of the evident relations, 



Z iZ 

^m,2 ^ —?= = Z/Q > 0, Vrn,l ^ —?= = ^t'O, 



the leading term of the expansion of DP in the small parameter, z/o/n <C 
1, can be expressed as 

-(2z/o + l)^(-z/o,0;4z/o)]+ (. 



Here ^E'(. . .) denotes irregular degenerate hypergeo metric function (cf 
(111)), whereas r(. . .) stands for the F-function 110 . 



Is — level 


2s — level 


uji = 3/16 


UJ2 = 1/8 


uji = 3/16 


UJ2 = 1/8 








1^12 


Z^21 


^22 


1^21 


V22 


^8/5 


^8/11 


2/v^ 


2/V5 




^8/5 


+ioo 


V2 



Table 1: The values of parameters Vij, i,j = 1,2. 

Let us apply the results obtained to particular photon energies adopted 
in the above-mentioned 1S-2S experiment in muonium atom [Z = 1). The 
latter is carried in the presence of two laser beams with the energies ui = 
3/16 a. u. (Ai = 244 nm) and uji = 1/8 a.u. (A2 = 366 nm). The corresponding 
values of Uij, i,j = 1,2 (|42|), ( ^31) and afs(co') are compiled in the Tables |l],^ 
In obtaining a£,(c<j2) we used the value. 



^(2s||r||^2) 



155.799140 - 46.045022i, 
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"fs(^) 








Ll)2 


-5.714105 


-4.962372 


29.853542 - 12.823175i 


89.818540 - 46.045022i 



Table 2: The values of 

being equal to the leading (77-independent) term in Eq.(^). Its imaginary 
part coincides, as it should, with 9q;£,(/2s) of Eq. (|5^) , whereas the real part 
defines the level shift at the photoionization threshold. According to Eq.(^), 
the numbers displayed enable to obtain, e.g. the following important dimen- 
sionless ratio: 

Here I^-^ = Pu^/Si, I^^ = P0J2/S2 are the beam intensities, Puii,oj2y 'S'1,2 being 
the corresponding powers and cross sections. The absolute value of the light- 
shift of the ns- level due to (one) photon beam of the field strength, E^^, is 
defined as (henceforth in this section the ordinary units are used) 

AeM = ]El^al{u), E^ ^ = 5.338- 10^^/^Eo mm-ns-^/^-mJ-^/^. 

4 \ ceo V S^T^ 

Here and stand for the energy of the beam within one pulse (in mJ) 
and the pulse duration (in ns); the beam cross section, S*^, is measured 
in mm?] Eq = m'le^/h'^ = 5.142 ■ 10^^ Volts/m denotes the atomic unit of 
electric field strength. For the typical values of these parameters adopted 
in the experiment {S^j = 2x3 mm^, = 28 ns, = 6 mJ) one obtains: 
Ei^ ~ 1.0 ■ 10~^£'o, i.e. the electric field employed happens to be rather 
weak. Accordingly, the average intensity within a pulse equals: to I^^ = 
Ap/{S^T^) ~ 3.57 ■ 10^ W/cm^. In the presence of two counterpropagating 
beams having the intensities, J^, J^, and the same frequency, u, the shift 

^The presence of two counterpropagating beams of the same frequency enables to avoid 
both the Doppler-broadening and the Doppler-shift of the Hne [Q . 
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of the level (at the frequency u) can be expressed in the form: 



= 4.6875 ( + I 3f^afs(^) ^ mm^ ■ ns"^ ■ mJ"^ ■ MHz. 



Here it is taken into account that a^g{uj), whose values are dislayed in the 
Table (|^), is measured in the units of ajj, ao = /{niee^) = 0.529 ■ 10~^ cm 
being the Bohr radius. As a result, the total energy shift between ns- and 
ms-levels caused by two counterpropagating beams of the same frequency uj 
is, then, given by 

( A'^ \ / 

= 4.6875 [-J^ + -JL-\ (^aliuj) - ^at{uj)) a^' mm' ■ ns"^ ■ 

An application of this formula to the case of 1S-2S experiment, assuming 
that parameters of the counterpropagating beams are identical {Suj-^ = S"^^ = 

= Si,^ = 2 X 3 mm^, 
7-^1 = = r^, = = 28 ns, = A;-^ = = A;-^ = 6 mJ), yields: 

A^aiK) = 11.9 MHz (59) 
AS2i{uj2) = i?^,^,A^2i(^i) =31.7 MHz. (60) 



3 Conclusion 

The numbers of Eqs.(p9D,(|60D may serve as a good illustration of the method 
employed. They are, however, of independent significance. The value of 
A£^2i(^^i) is in fair agreement with the result 0. One has to stress that in 
obtaining these energy shifts we used the average intensities. As it shown in 
[0], an account for a space-inhomogeneity of a laser field may increase each 
of these numbers by a factor of 10. It is relevant to point out here that the 
results obtained for = 1, 2 can be extended on the case of arbitrary n. Such 
type of calculation can be mostly efficiently performed on the basis of Eq.(|2^) 
valid for all n, e.g. with the help of Maple [l^, being an easy-to-use computer 
algebra program. Hence, in the contrast to a pure numerical scheme (say, 
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in a fashion of Beausoleil [^), one obtains analytical, rather than numerical, 
result which is already well adapted for further numerical calculation (one- 
fold integration). Such an integration is performed only at the very last stage 
of calculation and usually carried out in no time. Furthermore, by that means 
one considerably reduces numerical errors. Beyond that, analytical formulae 
admit straightforward computation of the various asimptotics with respect to 
all parameters encounted in them. As an example we can mention the result 
of Eqs.(|55D-(|56D when the photon energy, tu, tends to the threshold. Ins, being 



of particular importance for the problem under consideration. This argument 
may be considered as additional advantage of the method employed. Besides, 
as was already mentioned in the Introduction, the formula (^), as well as 
its particular cases, Eqs.(|5D|)-(^2|), enable a straightfoward calculation of the 
sums, Eq.(|]). In fact, S^^^uj) can be defined, in analogy with (^, as 

The analysis, whose details will be given elsewhere, shows that for low n it 
can be expressed in a closed form for any /i. The corresponding calculations 
prove to be simpler than those where the Green function is employed. We 
would interpret this circumstance, thereby, as a self-contained importance of 
the current method when it is applied to the problem under consideration. 
One would expect that this advantage will even be enhanced if the states 
with n > 4 are taken into consideration. It is interesting to emphasize 
here that one we can easily estimate the speed of the DP's decrease, as 
n — >• cx). Namely, according to Eg . (^) , a^(t^) x C/n^,n/i/o ^ 1- So that the 
threshold itself {E = 0), being the limit: Ens —^0, asn — oo, is not affected 
by the laser field. This is in agreement with the well known result due to Ritus 
[|T^ stating that the DP vanishes for continuum states. The ns— levels with 
n > 4 prove to be of particular importance owing to extensive experimental, 
as well as theoretical, investigation of the spectroscopic properties of the 
few-body systems, being carried out at present. In particuar, we consider 
the light-shift calculation in the 3s—, 4s—, 5s— levels as a subject of future 
publications. 
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5 Appendix 

A Orthogonality-relations 

Let us prove, say for n = 1,2, that the functions (pO|)-(|32|) satisfy "orthogo- 
nahty relation", Eq.([T3|). Namely, on calculating the integral. 



7'3/2 

{^PM = TTT^ I r'e-^'/^Mr)dr, n=l,2. 



we obtain: 



4 



After the variable changes, + 1) = t(z/ — l)/(z/ + 1), x/(a; + 1) = 

t{u - 2)/(z/ + 2), Eqs. ([OD , (|X^) take the final form: 

2^^VEu' /■! t'"" - 'uTi^) „ , 256 



Z3(zy-l)(z/+l)5(zy + 2)5io A (z.-l)(i.-2) ,\5 243Z3(z/2-l) 

1024 y6z/2 _ (2p||r||ls) 



243 Z3 (z/2 - 4) E- £2; 



(A.3) 



Here we have used the integrals: 



1 t^-ii-^.t 



(z.-i)(,.-2) ,\5 1296 z/4(z/-2) 



rft = ^- (A.6) 

2 \^ 4096z/4(z/- 2) ^ ^ 



5 

1/+2 



20 



These can be calculated by means of the following elementary relation: 
■ - ««) A = 1,±1^+1^±P±J1)(^ /' J^M. (A.7) 



{l-htf 4 6(1-6)4' 46 Jo (l-6t)4 

Eqs. ( [A.3| ),( [X^ ) imply that ifjiir), ip2{r) satisfy Eqs.(p!3|), as they should. 
The case of arbitrary n can be treated accordingly. 

B Matrix elements of i/jnir) 

Let us give here some details of a derivation of Eqs. (|3^ -([381). On calcu- 
lating the integrals, (ns| |r| n = 1,2, using Eqs.(pOD,(pT]) and explicit 
expressions 

we get: 

1, ,,,,,, fu -ly x^-''il + xY+'' , 2z/2 ^ 

-(ls||r||^i) = -32,^j^-^^(-^) I , ^ ._,5 (B.2) 



Z4(z/2 - 1)3 Vz/+ 17 Jo + Z^{y'^-\)' 
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After making the variable changes, + 1) = t(z/ — l)/(z/ + 1), + 1) 
t{v - 2)/(z/ + 2), Eqs.(P3),(P3) take the form: 



3 



1,„ 2iV /■! ti-'^fl-^t 

3 



^>-4)(. + 2)^yo 



^^^"^(^-^^^.t-i6-^(;^-^nB.5) 



2 
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Finally, on applying successfully recurrence relations, 



(1 - bty 



-dt 



62 Jo (l-6t)6 



dt + 2- 



-rlf-\- 

62 Jo [i-bty 

a2 /■! t^- 



,1-u 



dt 



1 



2-u-n /-i ti-'^ 



(l-6t)"+i n{l-b) 



n 



(l-6t)' 



+ 



-dt. 



62 Jo (1 - bty 



(B.7) 



together with ( IXTD to Eqs.(^),(gJ)' retrieve identities (|3|),(|3 
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Figure 2: The contour of integration in Eq.(32). 
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